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Grand Unification of the strong, weak and electromagnetic interactions into a single unified gauge 
group is an extremely appealing idea which has been vigorously pursued theoretically and exper- 
imentally for many years. The detection of proton or bound-neutron decays would represent its 
most direct experimental evidence. In this context, we studied the physics potentialities of very 
large underground Liquid Argon Time Projection Chambers (LAr TPC). We carried out a detailed 
simulation of signal efficiency and background sources, including atmospheric neutrinos and cosmo- 
genic backgrounds. We point out that a liquid Argon TPC, offering good granularity and energy 
resolution, low particle detection threshold, and excellent background discrimination, should yield 
very good signal over background ratios in many possible decay modes, allowing to reach partial 
lifetime sensitivities in the range of 10^* - 10^^ years with exposures up to 1000 ktonxyear, often 
in quasi-background-free conditions optimal for discoveries at the few events level, corresponding 
to atmospheric neutrino background rejections of the order of 10^. Multi-prong decay modes like 
e.g. p — > ij.~-K^ K'^ or p — > e"'"7r+7r~ and channels involving kaons like e.g. p K'^u, p — > e+A'" 
and p fj,'^ K° are particularly suitable, since liquid Argon imaging provides typically an order 
of magnitude improvement in efficiencies for similar or better background conditions compared to 
Water Cerenkov detectors. Up to a factor 2 improvement in efffciency is expected for modes like 
p e'^j and p — > fj.'^j thanks to the clean photon identification and separation from tt". Channels 
like p — > e^7r° or p ^ A'^'""") dominated by intrinsic nuclear effects, yield similar efficiencies and 
backgrounds as in Water Cerenkov detectors. Thanks to the self-shielding and SD-imaging prop- 
erties of the liquid Argon TPC, this result remains valid even at shallow depths where cosmogenic 
background sources are important. We consider the possibility of a very large area annular ac- 
tive muon veto shield in order to further suppress cosmogenic backgrounds at shallow depths. In 
conclusion, we find that this class of detectors does not necessarily require very deep underground 
laboratories, like those typically encountered in existing or planned sites, to perform very sensitive 
nucleon decay searches. In addition to a successful completion of the required R&D necessary to 
reach a relevant liquid Argon mass scale of 100 kton in a cost-effective way, we finally stress the 
importance of an experimental verification of the liquid Argon TPC physics potentialities to detect, 
reconstruct and classify events in the relevant GeV energy range. This experimental verification 
will require in addition to possible specific tests in charged particle beams, the collection of neu- 
trino event samples with statistics in the range of lOO'OOO events or more, accessible e.g. with a 
medium-sized detector at near sites of long baseline artificial neutrino beams. 

PACS numbers; 13.30.-a, 14.20.Dh, 29.40.Gx 



I. INTRODUCTION 

Grand Unification (GU) of the strong, weak and elec- 
tromagnetic interactions into a single unified gauge P, 0| 
is an extremely appealing solution which has been vig- 
orously pursued, theoretically and experimentally, for 
many years. An experimental hint in its favor is the 
apparent merging of the three coupling constants at a 
large energy scale (~ 10^^ GeV) when low energy mea- 
surements are extrapolated @. On the other hand, the 
most convincing experimental evidence for GU would be 
the direct observation of baryon number violation Q- 
The experimental search for decays of protons or bound- 
neutrons is therefore one of the most important and un- 
solved problem of particle physics. 

In the simplest Grand Unified Theories (GUT), nu- 
cleon decay proceeds via an exchange of a massive boson 
X between two quarks in a proton or in a bound neutron. 



In this reaction, one quark transforms into a Icpton and 
another into an anti-quark which binds with a spectator 
quark creating a meson. According to the experimental 
results from Super-Kamiokande 0, HI constraining the 
partial decay to t/B{p — > e"'"7rf) > 5.4 x lO'^^ years 
(90%C.L.), the minimal SU(5) 0, predicting a proton 
lifetime proportional to a~'^Mjr where a is the unified 
coupling constant and Mx the mass of the gauge boson 
X, seems definitely ruled out. 

Supersymmetry, motivated by the so-called "hierarchy 
problem" , postulates that for every SM particle, there is a 
corresponding "superpartner" with spin differing by 1/2 
unit from the SM particle In this case, the unification 
scale turns out higher, and pushes up the proton lifetime 
in the p — > c+tt*^ channel up to 10"^®^^ years, compatible 
with experimental results. At the same time, alterna- 
tive decay channels open up via dimension-five operator 
interactions with the exchange of heavy supersymmetric 
particles. In these models, transitions from one quark 
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family in the initial state to the same family in the fi- 
nal state are suppressed. Since the only second or third 
generation quark which is kinematically allowed is the 
strange quark, an anti-strange quark typically appears 
in the final state for these interactions. The anti-strange 
quark binds with a spectator quark to form a K meson 
in the final state [1]. The searches for decays p vK^ , 
n PK^, p — > /i+i^*^ and p — > e'^K'^ modes were also 
performed in Super-Kamiokande [lo| yielding counts 
compatible with background expectations, leadin g to lim- 
its on possible minimal SUSY SU(5) models [Til Il2l. [l3|. 
The theoretical predictions, however, vary widely, since 
there are many new unknown parameters introduced in 
these models. 

Other alternative models have been discussed in the lit- 

terature [H, M, M M M MMM,MMMM,M 

(see TableHl). In addition to the above mentioned GUTs, 
other supersymmetric SUSY-GUT, SUGRA unified mod- 
els, unification based on extra dimensions, and string- M- 
theory models are also possible (see Ref. [23 for a recent 
review). All these models predict nucleon instability at 
some level. Finally, it is also worth noting that theo- 
ries without low-energy super-symmetry 1231 f28| predict 
nucleon decay lifetimes in the range 10^^ years. 

Some experimental aspects of nucleon decay detection 
were discussed in Ref. [29]. Nucleon decay signals are 
characterized by (a) their topology and (b) their kine- 
matics. The presence of a lepton (electron, muon or neu- 
trino) in the final state is expected, and in general few 
other particles (two body decays are kinematically fa- 
vored), and no other energetic nucleon. The total energy 
of the event should be close to the nucleon mass and the 
total momentum should be balanced, with the exception 
of the smearing introduced by Fermi motion and other 
nuclear effects (nuclear potential, re-scattering, absorp- 
tion, etc.) for bound decaying nucleons. 

The search for nucleon decay therefore requires (1) ex- 
cellent tracking and calorimetric resolutions to constrain 
the final state kinematics and suppress atmospheric neu- 
trino backgrounds, (2) particle identification (in partic- 
ular kaon tagging) for branching mode identification (3) 
very massive detectors and (4) underground locations to 
shield against cosmic-ray induced backgrounds, although 
the exact required rock overburden depends on the cho- 
sen detection technology. Fine tracking in the low mo- 
mentum range (^100-1000 MeV/c) is fundamental for 
(dE/dx) measurement, particle identification and vertex 
reconstruction. 

In order to significantly improve current experimental 
results, next generation massive underground detectors 
satisfying the above requirements have to be considered. 
Given the variety of predicted decay modes open by the 
new theories, the ideal detectors should be as versatile as 
possible, very good in background rejection, and at the 
same time have the largest possible mass. The relevant 
factor is in fact Mxe, where M is the detector mass and e 
the signal detection efhciency after cuts to suppress back- 
grounds, which depends on the considered decay channel. 



Hence large masses must be coupled to fine tracking and 
excellent calorimetry, to suppress atmospheric neutrino 
and cosmogenic backgrounds with a good signal selec- 
tion efficiency. Furthermore, the detector should be sen- 
sitive to several different channels in order to better un- 
derstand the nucleon decay mechanism. Since there are 
about 6 X 10"^^ nucleons per kton of mass, the proton 
lifetime limit (90% CL) in case of absence of signal and 
backgrounds is about Tp/B > M (kton) x e x T xlO^^ 
years, where T is the exposure in years and B the as- 
sumed branching fraction for the searched mode. There- 
fore, the required effective mass Mxe to reach 10^^ years 
is in the range of 100 kton assuming T=10 years. 

Such massive underground detectors will be sort of ob- 
servatories for rare physics phenomena like astrophysi- 
cal neutrino detection and nucleon decay searches, with 
possible synergies with existing or new artificial neutrino 
beams for improved understanding of neutrino flavor os- 
cillations, including the possible identification of direct 
CP-violation in the leptonic sector (see e.g. Ref. [30| and 
references therein). 

Among the various options and technologies currently 
thought of (see Ref. ^ and Refs. |1 HI [H, Hi, 
IsgI . IstI. IssI ), the Liquid Argon Time Projection Cham- 
ber ^,140, IM El Ss, 44, 44111, 113, m (LAr TPC) is a 
powerful detector for uniform and high accuracy imaging 
of massive active volumes. It is based on the fact that in 
highly pure Argon, ionization tracks can be drifted over 
distances of the order of meters. Imaging is provided 
by position-segmented electrodes at the end of the drift 
path, continuously recording the signals induced. Tq is 
provided by the prompt scintillation light. 

Early work on the detection of nucleon decays in liquid 
Argon can be found in Ref. [i^. The liquid Argon TPC is 
a very promising detector option which satisfies at best 
the above mentioned requisites in terms of granularity, 
energy resolution and imaging capabilities if it can be 
extrapolated to the relevant mass scale. 

In this paper we study the performance of a very mas- 
sive Liquid Argon TPC as a nucleon decay detector. In 
the present study, we address for the first time the ef- 
fect of the charged cosmic rays background as a function 
of the depth (i.e. rock overburden) of the underground 
detectors. In particular, we address the possibility to 
perform proton decay searches in "shallow depth" con- 
figurations [si'l. Since the new large underground detec- 
tors will require either a new site to be excavated or the 
extension of an existing infrastructure, it is important 
to understand if very-sensitive nucleon decay searches 
do necessarily require deep underground locations, like 
those ty pica lly encountered in existing or planned labo- 
ratories |50 |. 

The paper is organized as follows: some considerations 
about the conceptual design of the detector are outlined 
in sectionini The discussion of the assumed detector con- 
figuration and the simulated physics process for signal 
and background events are described in section [TTTl The 
description of the analysis cuts designed to suppress at- 
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Model 


Ref. 


Modes 


TN (years) 


Minimal SU{5) 


Georgi, Glashow 


P 


e+n" 




Minimal SUSY SU{5) 


Dimopoulos, Georgi fill. Sakai [12] 


P 


VK+ 






Lifetime Calculations: Hisano, 


n 


vK° 


1028 _ 1032 




Murayama, Yanagida [13] 








SUGRA 517(5) 


Nath, Arnowitt [14, JJ,] 


P 


vK+ 


1032 _ 103-1 


SUSY SO{W) 


Shafi, Tavartkiladze [16] 


P 


uK+ 




with anomalous 




n 


uK° 


1032 _ 1035 


flavor [/(I) 




P 






SUSY SO{W) 


Lucas, Raby [17], Rati [18] 


P 


vK+ 


1033 „ 1034 


MSSM (std. d = 5) 




n 


uK° 


1032 _ 1033 


SUSY 50(10) 


Pati [18] 


P 


vK+ 


1033 _ 1034 


ESSM (std. d = 5) 








< iqSs 


SUSY S'O(10)/G(224) 


Babu, Pati, Wilczek [19, 20, 21], 


P 


vK+ 


< 2 ■ 10^^ 


MSSM or ESSM 


Pati fl8] 


P ~* 






(new d = 5) 






B ' 


- (1 - 50)% 


SUSY SU{5) or SO{W) 


Pati fl8] 


P 






MSSM (d = 6) 










Flipped SU{5) in CMSSM 


Ellis, Nanopoulos and Wlaker[22] 


P 




^q35 _ ^q36 


Split 517(5) SUSY 


Arkani-Hamed, et. al. [23] 


P 




1035 _ 1037 


517(5) in 5 dimensions 


Hebecker, March-Russell[24l 


P 
P 




1034 _ 1035 


SU{5) in 5 dimensions 


Alciati et.al. [25] 


P 




1036 _ 1039 


option II 










GUT-like models from 


Klebanov, Witten[26] 


P 


e+n" 


~ 10^^ 


Type IIA string with D6-branes 











TABLE I: Summary of the expected nucleon lifetime in different theoretical models. 



mospheric neutrino and muon-induced backgrounds from 
proton and neutron decay channels are given in sec- 
tion |TVl Finally, the results obtained in terms of lifetime 
sensitivities for the considered decay modes are given in 
section |Vl 



II. DETECTOR CONCEPT 

Our analysis assumes the scalable concept of a liquid 
Argon TPC, as proposed in [Flj (see Figure [1]). Other 
designs have been presented in Ref. [11]. An LOI based 
on a more standard configuration and a surface detector 
has also been submitted to FNAL [ssj . 

The design of Ref. [5l| relies on (a) industrial tankers 
developed by the petrochemical industry (no R&D re- 
quired, readily available, safe) and their extrapolation 
to underground or shallow depth LAr storage, (b) novel 
readout methods with embedded charge chain for very 
long drift paths, based on e.g. LEM readout, (c) new 
solutions for very high drift voltage, (d) a modularity 
at the level of 100 kton and (e) the possibility to em- 
bed the LAr in a magnetic field [13] ■ Such a scalable, 
single LAr tanker design is the most attractive solution 
from the point of view of physics, detector construction. 



operation and cryogenics, and finally cost. 

An R&D program is underway with the aim of opti- 
mizing the design ^55j . Concerning the further consoli- 
dation of the technology, we point out that a test-beam 
dedicated to the reconstruction and separation of elec- 
trons from neutral pions has been discussed [56| . In ad- 
dition, a ~ 100 ton liquid Argon TPC to complement the 
1 kton Water Cerenkov detector at the potential 2 km 
site 2.5° off-axis from the T2K beam has also been con- 
sidered [srj . If realized, this unique experimental setup 
will allow to compare the performance of the liquid Argon 
TPC to the Water Cerenkov ring imaging and to recon- 
struct neutrino events directly in the same beam with a 
statistics of more than 100 '000 events per year, sufficient 
to extrapolate the atmospheric neutrino background in a 
potential 1000 kton x year exposure. ICARUS T600 
to be commissioned in the coming years, will detect too 
few contained events to demonstrate the ultimate physics 
performance of the technology. 

When operating TPCs with long drift paths and charge 
amplification at the anode as the one considered here, 
one must pay attention to possible drift field distor- 
tions caused by the space charge created by the posi- 
tive ions slowly drifting towards the cathode. Indeed, 
it is well known that ion (and hole) mobility in liquid 
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FIG. 1: Conceptual design of a lOOkton LAr TPC detector (see Ref. [Sjl for details). To suppress cosmogenic backgrounds, 
we consider the possibility to locate three large area veto planes to detect remaining penetrating cosmic ray muons passing in 
the vicinity of the Argon detector (see text). 



Argon is small compared to that of electrons {fiions ~ 
(10^^ 10^^)/Xe), and inconsistencies among experimen- 
tal measurements were attributed to the influence of the 
motion of the liquid itself [Ell . 

The problem of ions accumulating in the TPC drift 
volume has been effectively overcome by introducing an 
additional grid (^ate') between the drift volume and the 
readout plane [53, [S^l • The gate is normally closed and 
only opens for a short period of time when an external 
trigger signals an interesting event. This method is easily 
applicable to our design in the case of a pulsed artificial 
neutrino beam source, as the one considered in Ref. [soj . 

In the case of proton decay (or atmospheric neutrino 
background) where the detector must be continuously 
sensititive, the prompt scintillation light with an appro- 
priate detection threshold (e.g. > 100 MeV) could be 
used to signal the (rare) occurrence of interesting events. 
Unfortunately, this is impracticable at shallow depths be- 
cause, as it will be shown later, the fiducial volume oc- 
cupancy is dominated by remaining cosmic ray events 
crossing the detector: at a rock overburden < 1 km w.e., 
the average number of cosmic ray muons crossing the 
detector within a drift time is > 1 (see Table IIIip . This 
situation worsens when the rock overburden is decreased. 

Hence, either the rock-overbuden is > 1 km w.e. or the 



TPC must be continuously active and we must consider 
space charge effects from accumulating ions in the liquid 
Argon volume. It is possible to show that asymptotically 
the stationary density of ions pi in the drift volume of a 
double phase readout system with charge amplification 
is given by: 



Pi{t oo) 



(1) 



where G is the average gain of the charge readout system, 
£ feedback is the efficiency for ions to feed-back into the 
gas region below the charge amplification planes, ri^-i is 
the efficiency for ions to be transmitted from the gaseous 
phase into the liquid phase, egrid is the transparency of 
the extraction grid when traveling from the higher ex- 
traction field to the lower drift field, < dE/dx > is the 
average stopping power of the cosmic ray muons, x is the 
average path of the muons in the fiducial volume, (/i^ the 
flux of cosmic muons at the detector upper surface per 
unit surface and per unit time, fj-ions is the ion mobility, 
E the electric drift field and We{E) the field-dependent 
average energy expended by a cosmic ray muon to cre- 
ate an ion-electron pair in the medium (field dependence 
arises from columnar ion-electron recombination [46j). 
Assuming a stationary situation, ionization electrons 
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will drift in a medium homogeneously filled with the 
above calculated density of ions. The resulting free elec- 
tron lifetime can be estimated as Tg = 1/ {kr{E)pi) where 
kr is the electron-ion recombination rate (which depends 
on the electric field). 

We note that e feedback ^ 1, because of electron diffu- 
sion in the amplification gap, as a result of which ions 
can follow field lines ending up on the electrodes of the 
amplification system rather than on the external field 
lines; this is particularly true for the devices with multi- 
ple stages that we operate [6l|. In addition, e^rid ^ 0.3 
since the extraction field is typ. 3 kV/cm and the drift 
field is typ. 1 kV/cm. 

The efficiency rjy-i for the transfer of positive ions from 
the vapour into the liquid Argon phase has, as far as we 
know, not being measured. When the distance between 
the ion and the liquid- vapor interface is greater than sev- 
eral Angstroms (the Argon bond length is 3.7 A), the 
liquid may be treated as a continuum (6^ . Assuming a 
planar interface between the liquid and vapor, the situ- 
ation is then comparable to a boundary-value problem 
with dielectrics and can be solved with the charge-image 
method in a single dimension: a point charge q in the 
medium e at the interface with a medium e' feels a charge 
image q' in the medium e' given by q' = — (e' — e)/(e'+e)g. 
For e' > e, q.q' < and the force is attractive. For e' < e, 
q.q' > and the force is repulsive. Given the difference 
in dielectric constants between the liquid Argon e/ ~ 1.54 
and its gas phase ~ 1, a charge in the liquid will feel 
a repulsive force from the interface, while a charge in 
the gas will be attracted to the interface and will tend to 
"stick" at the surface. The potential energy of the charge 
q placed at a distance d > from the planar liquid-vapor 
interface is [11 [fil [fil [fii : 
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where the index corresponds to the charge in the liquid 
{I) or in the vapor (v). Classically, the potential barrier is 
infinite at the interface {d = 0) and a quantum mechani- 
cal treatment is needed. In practice, one can also assume 
that in order to cross the boundary the charge must pass 
above a finite, classical potential, whose height is to be 
determined experimentally [6^. 

In the presence of an external difference of potential to 
create a drift field, the total potential energy is V{x) = 
Vi{\x\)+qEix iov X < (liquid phase) a,ndV{x) = Vy{x) + 
qEyX for X > (vapor phase), where Ei (resp. E^ = 
{ei/ey)Ei) is the induced electric field in the liquid (resp. 
vapor) phase. 

Experimental data has been collected on the extrac- 
tion of quasi-free electrons from liquid into gaseous Ar- 
gon [g^. The potential felt by such electrons is illus- 
trated in Figure [21 It is plotted in units of kT where 
T = 87 K is the temperature of the liquid Argon at nor- 
mal pressure. Photoelectric effect measurements indicate 
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FIG. 2: Illustration of the potential felt by quasi-free electrons 
as a function of the position a:: in A across the liquid-vapor 
interface in the presence of an external drift field. The poten- 
tial is given in units of kT where T = 87 K is the temperature 
of the liquid Argon at normal pressure. 
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FIG. 3: Same as Figure [2] but for positive ions. 



that the electrons in the liquid are shifted by an amount 
Vb = —0.21 eV with respect to the energy in vacuum (or 
vapor) [66j . The external drift field and the repulsive po- 
tential at the interface generate a minimum potential in 
the liquid phase at a distance d'^ = {Aij eEi)^!"^ from the 
interface with a value V^^ = 2(y4;ei?i)^/^ and similarly in 
the gas phase. The barrier height as seen from the min- 
imum in the liquid at x = —dm is therefore AT^(i?) = 

v.-vl-v^^ ^ v,^ 2eV2 [a]I^e\i^^a]J^eII-'). 

The mean life time to traverse the interface (or trapping 
time) can be analyzed with the aid of the Smoluchowski 
equation [63j . predicting Te oc exp(AI^/fcT). Experi- 



mental data however indicate that the above-calculated 
thermionic current is increased by a rate oc ^E/Xi where 
Ai is the momentum-transfer mean free path, as pre- 
dicted by the Shottky model of electric field enhanced 
thermionic emission (65j . giving t^, oc exp(AV^/fcr) /E, 
with a stronger influence of the drift field. 

In the case of the ions, the potential for a positive 
charge is illustrated in Figure[31 The ion is attracted very 
near the interface and must tunnel through the repulsive 
potential existing on the liquid side of the interface in or- 
der to enter the liquid. A semi-classical approach is not 
justified. The transport of Tl positive ions from xenon 
vapor to liquid was successfully reported with a ~ 10% 
efficiency in Ref. [131 ■ There appears to be a difference 
in the behavior of positive ions between the liquid-to-gas 
and gas-to-liquid transition While electrons could 

traverse the liquid- vapor interface the positive ions could 
not. This is consistent with the fact that the ion mobility 
in the liquid is too small to allow for the external field to 
enhance the thermionic emission. However, the transmis- 
sion of negative ions through the liquid- vapor interface of 
neon has been observed with trapping time of the order of 
10-100 s at fields of 0.1-1 kV/cm [ej]. This is explained 
by the type of clustering around the different species of 
particles. If a sort of bubble surrounds negative charges 
then the breaking of the bubble at the surface could let 
them escape. In conclusion, the phenomenon of transfer 
of the ions through the interface is expected to be rather 
complicated, not well understood and has not been mea- 
sured. In the case relevant to this paper, it is likely that 
some fraction of the Argon ions will cross into the liquid, 
although possibly with long trapping times, and experi- 
mental studies are needed to assess how much. For safety, 
we will conservatively assume full transmission of Argon 
ions across the interface. 

We note that competing with the process of transfer 
from the gas to liquid is the radial drifting of ions at the 
surface of the liquid-gas interface towards the edges of 
the volume, until they reach the field shaping electrodes. 
There they are neutralized by the external power supply. 

We now compute the free electron lifetime due to the 
recombination with the accumulating ions. In the worst 
case considered in this paper corresponding to the "un- 
der the hill location" (see Table IIIip . we contemplate 
~ 1000 muons crossing the detector per second. We con- 
servatively assume that each of these muons will verti- 
cally cross the entire drift region, hence a: ~ /i = 20 m, 
that the average charge gain is G ~ 150 (see Ref. [5l|). 
very conservatively that tfeedhack — 0.5, -q^-i ~ 1 al- 
though we expect this to be a pessimistic assumption, 
Cgrid — 0.3 and that the ion mobility in the range 
(0.2 1) X 10"^ cm^/V/s (53]. These assumptions yield 
Pi ~ (0.7 3.7) X 10^ ions/cm'^. Assuming a recombi- 
nation rate = 1 x 10~^ cm^/s, consistent with mea- 
surements at the relevant electric drift fields , implies 
that the electron lifetime due to ion accumulation in the 
medium is Tg > 30-^ 140 ms, to be compared with a max- 
imum electron drift time ~ 10 ms. We stress that actual 
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FIG. 4: Drift field map with an unrealistic ion charge density 
of 2 X lO'^ ions/cm^=4 x 10"* C/m^ (see text). The field 
distortion is ±30%. Dark blue corresponds to 0.5 kV/cm and 
red to 1.5 kV/cm. 
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FIG. 5: Same as Figure U with a ion charge density of 2 x 
10" ions/cm^ =4 x lO"" C/m^ The field distortion is ±3%. 



values of e feedback and r]jj„i are expected to increase these 
values even further. We conclude that in all configura- 
tions considered in this paper, charge attenuation due to 
recombination with ions in the medium is expected to be 
negligible compared to other charge attenuating effects, 
like e.g. attachment to electronegative impurities in the 
Argon [45]. 

We now turn to the drift field map distortions. The 
results of numerical simulations for the electric field as- 
suming a homogeneous positive charge density distribu- 
tion in the liquid Argon fiducial volume of resp. pi = 
2 x 10^ ions/cm"^ and 2 x 10* ions/cm"^ and the contribu- 
tions of ions directly produced in the ionization process 
are shown in resp. Figures |4] and [3 The cathode is placed 
at —2 MV and the volume is enclosed by field shaping 
electrodes. Dark blue corresponds to 0.5 kV/cm and red 
to 1.5 kV/cm. With an unrealistic ion charge density of 
2 X 10^ ions/cm^=4 x 10~^ C/m^, the field is distorted by 
±30%. As expected, for a more realistic charge density 
of 2 X 10* ions/cm^, field variations are ±3%. The effect 
of ions directly produced in the ionization is negligible. 

We point out in this context that the realization and 
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successful operation on surface of a 5 m long detector 
column ^55,] (ARGONTUBE) will allow to experimen- 
tally verify these hypotheses and to prove the feasibility 
of detectors with long drift paths, hence representing a 
very important milestone in the conceptual proof of the 
detector design. 

III. SIMULATION FRAMEWORK 

As already pointed out in Section |T1 the nucleon decay 
searches aim at discoveries at the few event level in an 
atmospheric neutrino background sample of more than 
10^ events. The understanding of this background will 
require the collection of a commensurate sample in a neu- 
trino beam of the relevant energy. 

In the meantime, a reliable MC description of all phys- 
ical effects involved is important to correctly estimate 
nucleon decay sensitivities. Whenever possible, our re- 
sults are based on full simulations. Several neutron and 
proton decay channels with proper treatment of nuclear 
effects have been simulated. Likewise, the two dominant 
background sources relevant to our search were simulated 
in details: the atmospheric neutrino background and the 
cosmic muon-induccd background. The next sections de- 
scribe this in more details. 



A. Detector geometry and simulation 

We considered as reference detector a homogeneous 
volume of liquid Argon in a single giant volume composed 
of a cylindrical volume of 20 m height and 70 m diameter 
for a total of 100 kton. This geometry was implemented 
in a full simulation based on the GEANT4 toolkit [MfTOt. 

For each event generated within the liquid Argon vol- 
ume, final state particles are transported through the 
medium, with the possibility of secondary interactions. 
The detector effects have been included in the produc- 
tion and transport of the events by simulating the liquid 
Argon response. For ionizing particles, the deposited en- 
ergies in microscopic sub-mm'^ volumes of Argon are used 
to compute the amount of ionization 46] and scintillation 
yields |7i[. Ionization charge is collected in two perpen- 
dicularly segmented striped readout planes with 3 mm 
pitch. This information is digitized taking into account 
the response of a typical potential charge readout elec- 
tronic preamplifier and is used to create the two perpen- 
dicular charge readout views. Scintillation light is propa- 
gated through the medium and collected on photodetec- 
tors located on the inner surface of the detector 1961. 



B. Nucleon decay signal simulation 

Several neutron and proton decay channels have been 
studied (see Section IIV[) . For each channel, 2000 sig- 
nal events were generated and fully simulated inside the 



detector. Nucleon decay events are characterized by a 
definite value of the total energy and by the fact that 
the total momentum of the decay products must be zero. 
These features, which are true for the decay of a free 
nucleon, are only approximately verified for a nucleon 
bound in a nucleus. 

The nuclear effects, the distortions of the energy and 
momentum distributions due to Fermi motion (since the 
recoil nucleus is not measured), and the reinteraction of 
decay particles with the nucleus have been treated with 
the FLUKA 72] package. The treatment is similar to that 
used for nuclear effects in neutrino interactions, see e.g. 
Ref. j/zj], and will be used also in the the context of 
simulation of atmospheric neutrino background (see Sec- 
tion UTTC]). 

In a given signal event, energy and relative orientation 
of the decay products are selected according to conserva- 
tion laws and phase space. The resulting configuration is 
inserted into the target nucleus and used as initial step 
for the PEANUT [tI] simulations. The position of the de- 
caying nucleon in the nucleus is sampled from a probabil- 
ity distribution proportional to the density profile for the 
selected nucleon type. The Fermi motion of the original 
nucleon is then sampled according to the local Fermi dis- 
tribution. The decay products are followed like any other 
secondary particle in PEANUT, thus exploiting all the de- 
tails of its nuclear modeling, as described in the following 
paragraphs. See Ref. 72] and references therein for full 
details. 

Nuclear effects in nucleon decays can be roughly di- 
vided into those of the nuclear potential and those due 
to reinteractions of decay products. Bound nucleons in 
nuclei are subject to a nuclear potential. The Fermi en- 
ergy (or momentum) must be calculated from the bottom 
of this nuclear potential well, and the removal of a nu- 
cleon from any stable nucleus is always an endothermic 
reaction. When a nucleon decays, some energy is spent 
to take it out of this well: the minimum energy is given 
by the nucleon separation energy (around 8 MeV), and 
corresponds to the decay of a nucleon at the Fermi sur- 
face. In this case, the daughter nucleus is left on its 
ground state. More deeply bound nucleons decay leaving 
a hole in the Fermi sea, that corresponds to an excita- 
tion energy of the daughter nucleus, and an additional 
loss of energy of the decay products. This energy is then 
spent in evaporation and/or gamma deexcitation. Thus, 
the invariant mass of bound nucleon decay products is 
expected to be always slightly smaller than the free nu- 
cleon mass, and spread over a range of about 40 MeV. 
Correspondingly, the Fermi momentum of the decaying 
nucleon is transferred to the decay products and compen- 
sated by the recoil of the daughter nucleus. Additional 
momentum distortions come from the curvature of par- 
ticle trajectories in the nuclear potential. 

Reinteractions in the nuclear medium also play an im- 
portant role. Decay products can lose part of their en- 
ergy in collisions, or even be absorbed in the same nu- 
cleus where they have been created. This is particularly 
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true for pions, that have an important absorption cross 
section on nucleon pairs, while kaons have smaller inter- 
action probability. 

Nucleon-nucleon total cross sections, both elastic and 
inelastic, used in FLUKA are taken from available exper- 
imental data. Elastic scattering is explicitly performed 
according to the experimental differential cross sections. 
Pion induced reactions are more complex, mainly because 
of two- and three-nucleon absorption processes. Above 
the pion production threshold, the inelastic interactions 
are handled by the resonance model. Other pion-nucleon 
interactions proceed through the non-resonant channel 
and the p-wave channel with the formation of a A res- 
onance, with modified resonance parameters taking into 
account nuclear effects. 

The conservation of strangeness leads to very different 
interactions of the K+, K° and K^, K" with nucleons 
at low energies. K~ have a large cross section for hy- 
peron production, with the Stt and Att channels always 
open. A detailed treatment of interactions has been 
developed in PEANUT, whereby the nucleon system 
is described by phase shift analysis. This treatment was 
used in the present analysis, where it was very relevant 
for the nucleon decays involving kaons. In particular, 
involved in proton decays are below 700 MeV/c, where 
the inelastic contributions are very small and elastic scat- 
tering does not affect strongly the efficiency Charge ex- 
change +n ^ K'^ + p can lead to a loss of efficiency. 
In the simulations of p — > vK^ decays, less than 4% of 
the kaons were lost due to nuclear processes (see Sec- 
tion |IV|. In comparison, in the simulated p — > e+7r'^ 
decays, the neutral pion is absorbed within the nucleus 
with a probability of ~ 50 %. 

In the treatment of nuclear reinteractions, the cross- 
sections are modified to avoid too short mean free paths 
in nuclear matter by taking into account Pauli blocking, 
antisymmctrization of fermions, hard-core effects and for- 
mation zone or coherence length for non- fragmenting pro- 
cesses. The intra-nuclear cascade step goes on until all 
nucleons are below 50 MeV and is followed by a pre- 
equilibrium emission phase in which additional nucleons 
can be emitted. 

Many benchmarks showing good agreement between 
PEANUT results and nuclear reaction data can be found in 
the main Ref. (t^ . 



C. Atmospheric neutrino background 

The atmospheric neutrino flux has been computed by 
several groups. We take results from the FLUKA group [Tj 
and the HKKM 2004 ^ (Honda) models. The FLUKA 
model is estimated to have a 7% uncertainty for the pri- 
mary spectrum, 15% for the interaction model, 1% for the 
atmosphere profile, 2% for the geomagnetic field and a 
total 17% uncertainty The HKKM model has about 10% 
uncertainty when the neutrino energy is below 10 GeV, 
but the uncertainties are still large when the neutrino 



energy is above 10 GeV, due to the uncertainties of the 
primary cosmic ray flux and the interaction model above 
100 GeV. 

The neutrino-cross sections and the generation of neu- 
trino interactions are based on the NUX [t^ code. This 
code can be used in fixed target mode for incoming neu- 
trino energies from ~ 20 MeV up to 10 TeV. NUX can 
generate neutral current (NC) and charged current (CC) 
processes, and includes all basic processes to properly 
describe neutrino interactions in different kinematical 
regimes. In particular, a distinction is made among the 
following processes: quasi-elastic process (QE), inelastic 
process (DIS) [131 and charm production. These pro- 
cesses are appropriately matched to reproduce the total 
inclusive and existing exclusive experimental data to the 
best of current experimental knowledge. The mass of the 
outgoing lepton and of the target nucleon are taken into 
account in the kinematics and in the differential cross- 
section. In order to take into account nuclear effects, an 
interface with FLUKA was implemented, yielding the so- 
called NUX-FLUKA generator: a primary nucleon is cho- 
sen according to the density profile for the selected nu- 
cleon type and the Fermi motion of the original nucleon 
is sampled according to the local Fermi distribution, just 
as in the case of the nucleon decay simulation (see Sec- 
tion |TlIB|. A neutrino interaction is generated and the 
final state intra-nuclear cascade is taken into account by 
propagating final state particles through the nucleus with 
PEANUT. 

The NUX-FLUKA model was benchmarked with the data 
of the NOMAD [t^I experiment at high neutrino ener- 
gies (> 5 GeV), and was compared at low neutrino en- 
ergies (< 5 GeV) with other available generators like 
NUANCE [tH] and was shown to give very similar par- 
ticle multiplicities and distributions (79j . The produc- 
tion of strange particle at high energy was studied in 
NOMAD (s^l and found to be properly described after 
tuning of the JETSET 'si'l fragmentation parameters. At 
low energy, the production of kaons primarily proceeds 
through the decay chain of baryon resonances. In this 
regime, the yield of kaon is very strongly suppressed com- 
pared to that of pions. 

Table HIl reports the computed number of atmospheric 
neutrino events per neutrino fiavor and per process type, 
normalized to an exposure of 1 ktonxyear. In or- 
der to reach nucleon decay sensitivities in the range of 
2^g35 years, exposures in the range of 1000 ktonxyear 
will be considered. Therefore, the simulated data sets 
correspond to a total exposure of 1000 ktonxyear or 
about 250'000 background events. Part of this atmo- 
spheric neutrino background sample was fully simulated 
and fully reconstructed. Typical atmospheric neutrino 
events in Argon viewed by two-dimensional charge read- 
out planes with a pitch of 3 mm are shown in Figures [H] 
and[7l Full reconstruction of atmospheric neutrino events 
at low multiplicity is rather straight-forward. In partic- 
ular, quasi-elastic events reconstruction is discussed in 
Ref. [82l | . Actual performance on real quasi-elastic events 
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FIG. 6: Typical i/^ and i/g QE event in liquid Argon detector {i/^^ + n ^ p + fi and i^^ + n ^ p + e 




FIG. 7; Typical Vfj, NC event in liquid Argon detector (i/^ + p ^ V/j. + p + tt") 



data can be found in Ref . [43| . See also Ref. [83| for a dis- 
cussion of precision atmospheric neutrino measurements 
with a liquid Argon TPC. 

In order to have the largest sample of events to es- 
timate the backgrounds after selection cuts, a fast re- 
construction program was developed which used Monte- 
Carlo tracks as input for pattern recognition. Energy and 
angular resolution of the detector were properly taken 
into account by means of the full simulation. On the 
other hand, particle reconstruction efficiency was intro- 
duced as momentum thresholds which depended on the 
particle type. These latter were approximately 5 MeV/c 
for electrons, 10 MeV/c for photons, 20 MeV/c for muons 
and pions, 30 MeV/c for kaons, and 300 MeV/c for pro- 



tons. Checks were performed to ensure that the fast 
simulation reproduced the general features of the back- 
ground events on the subsample of fully reconstructed 
events. 



D. Cosmic muon-induced background 

A second relevant source of background comes from 
particles produced by cosmic muons crossing and inter- 
acting in the liquid Argon volume, or muons avoiding the 
active liquid Argon but interacting in the vicinity of the 
detector. The main physical parameter determining the 
rate of this background is the rock overburden. 
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Atmospheric 


Ve, CC (evts/kt/year) 


Ve. CC (evts/kt/year) 


Flux Model 


DIS QE Tot 


DIS QE Tot 


FLUKA 


22.9 49.0 71.9 


7.2 8.4 15.6 


(B(GeV)) 


(2.83) (0.62) 


(3.12) (0.79) 


HKKM 


22.7 45.6 68.3 


7.2 7.8 15.0 


(£(GeV)) 


(3.01) (0.64) 


(3.42) (0.85) 


Atmospheric 


CC (evts/kt/year) 


CC (evts/kt/year) 


Flux Model 


DIS QE Tot 


DIS QE Tot 


FLUKA 


42.0 80.1 122.1 


15.3 18.5 33.8 


(£(GeV)) 


(4.53) (0.72) 


(4.89) (0.90) 


HKKM 


43.3 74.0 117.3 


15.8 17.4 33.2 


(S(GeV)) 


(4.78) (0.77) 


(5.19) (0.98) 


Atmospheric 


V NC (inelastic) 


V NC (inelastic) 


Flux Model 


(evts/kt/year) 


(evts/kt/year) 


FLUKA 


23.2 


9.0 


(S(GeV)) 


(3.73) 


(4.08) 


HKKM 


23.4 


9.1 


(£(GeV)) 


(3.97) 


(4.40) 



TABLE II: QE, DIS event rates per kton per year, and average 
energies for FLUKA 2002 and HKKM 2004 flux models and 
NUX neutrino cross-sections. In order to reach nucleon decay 
sensitivities in the range of 10'^^ years, exposure in the range 
of 1000 ktonx year will be considered. 



The cosmic muon intensity as a function of under- 
ground depth is obtained from Crouch's [13] fit to the 
world data and further updated by the Particle Data 
Group jssj l. This parameterization agrees to better than 
6% with the most recent MACRO measurements [s^ at 
Gran Sasso. However, it cannot be apphed to the case of 
a shallow depth detector. For this purpose, we generated 
cosmic ^"""'s and /i~'s according to double differential dis- 
tributions found in Ref. [ssj on a large hemisphere at the 
surface of the Earth and transported them through rock. 
In the simulation, all physical processes were taken in 
account, including multiple scattering and energy losses. 

We aimed at giving here the general trend and the 
magnitude of the background, while the precise deter- 
mination of the muon induced background will depend 
on many site-specific features: (a) rock chemical com- 
position, (b) exact distribution of rock overburden (c) 
detector configuration, cryostat, distance to rock, ... 

Many relevant processes depend on rock chemical com- 
position. For example, the neutron production per unit 
of rock mass depends on the average atomic number. 
The expected dependence according to Ref. [s^l is dN oc 
means a difference of about 25% in neu- 
tron production going from standard rock ((A) — 22) to 
the salt rock ((A) = 30). Particle propagation through 
rock of secondary particles produced in muon interac- 
tions is also influenced by chemical composition, in first 
approximation with the same behaviour as the nuclear 
interaction length. 



The angular distribution of muons underground de- 
pends on the actual distribution of rock overburden. If 
the Earth surface above the site is approximately flat, 
the muon intensity at an angle d is given by the vertical 
intensity calculated at a depth h! — /i/cos0 multiplied 
by a factor 1 / cos 9 from the production probability in 
the atmosphere [84|. If, however, the detector is located 
in a mountain or under a hill, the muon intensity is de- 
termined by the actual rock overburden along any given 
direction. The most important effect is a different to- 
tal muon intensity with respect to the flat surface case, 
and a secondary effect is the difference in veto efficiency 
(see Section IIIIE[) due to the different muon direction- 
ality. Indeed, the deeper the rock overburden, the more 
vertical the surviving muons. 

Since, as we will show in the following paragraphs, sec- 
ondary particles that enter the detector are produced in 
the first one or two meters of material around it, the 
actual material distribution/composition of the LAr con- 
tainer may affect the background and the veto conditions. 
This effect has been neglected at this stage and we have 
assumed rock composition around the fiducial liquid Ar- 
gon volume. 

We have considered different underground depths and 
two different geometries (see Tables IIIII and lIVp in or- 
der to explore possibilities at shallow depths: (a) under 
fiat ground cover (e.g. a green-field site [111 or a mine 
with vertical access) at resp. approximately 3, 2, 1, 0.5, 
0.13 km water equivalent, corresponding to resp. 1130 m, 
755 m, 377 m, 188 m and 50 m of standard rock; and 
(b) detector under a hill configuration at shallow depth 
(see Figure H]). A green- field site, or a shallow-depth 
site under the hill were considered in the context of a 
synergy with an upgraded CNGS accelerator neutrino 
physics program [30'| . 

The simulations of muon-induced background were 
performed using FLUKA and GEANT4 toolkit and the re- 
sults were compared. They were carried out in three 
steps: (1) in the first step we estimated the number of re- 
maining cosmic muons crossing the detector or its vicinity 
per given unit time as a function of the detector's depth 
or its geographical configuration; (2) then we obtained 
particle fluxes at the detector surface from interactions 
of those cosmic ray muons in the surrounding rock as a 
function of the detector's depth; (3) the resulting neu- 
tron, kaon and lambda differential flux spectra, rescaled 
to the rate expected at a given depth, have then been 
used to simulate particle interactions inside the detector. 

As far as FLUKA simulations are concerned, the cross 
section of the photo-nuclear interaction of high-energy 
muons is taken from Ref. [HI and validation of the FLUKA 
models for muon transport and muon photo-nuclear in- 
teractions can be found in Refs. [87„ .89] . 

In the case of the detector under the hill, the topol- 
ogy of an existing location was studied [99[ and a simple 
but conservative description of the rock profile was im- 
plemented in a GEANT4 simulation. The parameters of 
the geometry in a 3D image from the GEANT4 simulation 
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Depth 

Water equiv. | Standard rock 


Code 


All 

Particles/s 


muons 

Particles/ 10 ms 


Particles/s 


> 1 GeV 
Particles/10 ms 


Effective mass 


Surface detector 


FLUKA 


1700000 


17000 


1300000 


13000 


- 


~ 0.13 km w.e. 


50 m 


FLUKA 


11000 


110 


10000 


100 


50 kton 


~ 0.5 km w.e. 


188 m 


FLUKA 


330 


3.3 


320 


3.2 


98 kton 




200 m 


GEANT4 


" 


" 


420 


4.2 


98 kton 


~ 1 km w.e. 


377 m 


FLUKA 


66 


0.66 


65 


0.65 


100 kton 


~ 2 km w.e. 


755 m 


FLUKA 


6.2 


0.062 


6.2 


0.062 


100 kton 


~ 3 km w.e. 


1.13 km 


FLUKA 


0.96 


0.01 


0.96 


0.01 


100 kton 


Under the hiU (see Figure [S| 


GEANT4 






960 


9.6 


96 kton 



TABLE III: Computed average number of muons entering the detector per unit time for various geographical configurations. 
The effective mass corresponds to the mass of Argon that can be used when, in both 2D readout views, a slice of size 10 cm 
around each crossing muon is vetoed. 



158 



200 m 




1 / 



70 m 
— ► 

5$ 20 



800 m 




FIG. 8: Under the hill detector configuration. The rock overburden at the vertical of the detector is 200 m and the shortest 
overburden is 158 m. This geometry has been implemented in a full simulation based on GEANT4. 



are shown in Figure [H The assumed rock overburden 
at the vertical of the detector is 200 m and the shortest 
overburden is 158 m. Cosmic ray muons were simulated 
on a circular area with a radius of 400 m at a height of 
200 m above the detector which corresponds to the top 
of the hiU. 

The first concern is the average number of muons en- 
tering the detector within a time interval of 10 ms, which 
corresponds to the assumed time for a full drift over 20 m 
at a field of IkV/cm ^Slj . The results are summarized in 



Table Uni 

The results from FLUKA and GEANT4 were compared 
considering a flat Earth surface profile and a detector at 
~ 200 m underground: muon rates agreed within a fac- 
tor of ~ 1.3. At depths deeper than about 1 km w.e. 
the rate of crossing muons is less than 1 per 10 ms. Shal- 
lower depths e.g. at rock overburdens less than 200 m are 
disfavored because of detector occupancy (at 50 m rock 
overburden the average rate of muons with more than 
1 GeV crossing in 10 ms is 100). If we assumed that a 
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slice of size 10 cm in a 2D readout view of the detector 
around each crossing muon cannot be used for physics, 
then the rate at resp. 50 m, 188 m would correspond to 
an available Argon mass of resp. 50 kton, 98 kton. We 
do not at this stage consider surface operation, as pro- 
posed in Ref. 53]. For the under the hill configuration, 
more than a million muons, corresponding to 156 drifts 
of 10 ms, were simulated and it was found that on aver- 
age 9.6 muons enter the Argon volume in 10 ms. This 
is a factor 2 higher than at an equivalent depth under 
flat surface, which is still very tolerable. 

We now focus on the second step where the rate and 
energy spectrum of remaining cosmic ray muons at the 
underground locations is used to study the number of 
particles entering the detector. We assume that events in 
which the parent muon enters (before or after the photo- 
nuclear interaction) the active LAr volume can be dis- 
carded thanks a veto based on the liquid Argon imag- 
ing. This leads us to restrict the background sources 
to neutral hadrons and, in particular, to neutrons^ neu- 
tral kaons and lambdas, produced either directly in muon 
photo-nuclear interactions or as secondary products in 
hadronic showers in the materials (e.g. rock) surround- 
ing the detector. We have so far conservatively neglected 
correlations between neutral hadrons penetrating inside 
the LAr in conjunction with (visible) charged hadrons. 
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FIG. 9: Energy distribution of neutrons entering the LAr 
volume as predicted by FLUKA. 



Neutral hadronic particles have been scored when en- 
tering the detector. The energy-integrated number of 
particles entering the detector per year of exposure at 
some depths are reported in Table IIVI Neutron, neutral 
kaon and lambda energy spectra on the detector surface 
are plotted in Figures [H [TO] and [TlJ Their shapes are 
practically independent of the depth of rock overburden. 
Since low energy neutrons are not a background source 
for nucleon decay searches, we will consider only neutrons 
with kinetic energy above 20 MeV. 

As expected, the particle yields per muons entering 
the detector, increase with depth, since the average muon 
energy increases as well. The neutrons yields in the FLUKA 
and GEANT4 simulations are consistent within a factor 
of 2. As far as neutral kaons are concerned, there is a 
difference of a factor of 10. In the case of lambda's, they 
were not found in the GEANT4 simulation while in FLUKA 
they are expected at the level of 0.04-0.05 times the rate 
of neutral kaons. In the following, we will use the yields 
from FLUKA to estimate backgrounds. 

In the third phase, neutrons, neutral kaons and 
lambda's were generated with their expected energy spec- 
trum and were transported in LAr until they underwent 
an inelastic interaction or decayed. These fully simu- 
lated interactions were subjected to the selection cuts 
and were used to estimate cosmogenic backgrounds (see 
Section HVC]). 
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FIG. 10: Energy distribution of neutral kaons entering the 
LAr volume as predicted by FLUKA. 



E. An active very large area cosmic muons veto 

The spatial distribution of photo-nuclear interactions 
in a 3 km w.e. depth from which at least one neutron 
enters the detector is plotted in Figure [TO] Most inter- 
actions occur along the detector vertical wall, at an av- 
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Configuration/ 
Depth 


Simulation 


Average number /I's 
entering LAr 
per 10 ms 


N< 

per year 


iutrons 
per fi in LAr 
per 10 ms 


Neuti 
per year 


"al kaons 
per /i in LAr 
per 10 ms 


Lambdas 


~ 0.5 km w.e. 

(188 m rock) 


FLUKA 


3.3 


1.9 xlO'' 


1.8x10"* 


4500 


4.3x10"'^ 


« 0.04 X Nko 


~ 1 km w.e. 
(377 m rock) 


FLUKA 


0.66 


5.5x10^ 


2.6x10-* 


1300 


6.2x10-'^ 


f» 0.05 X N^o 


~ 3 km w.e. 
(1.13 km rock) 


FLUKA 


0.01 


1.1 xlO"* 


3.6x10"* 


25 


8.2x10"^ 


f» 0.06 X N^o 


Under the hill 
(see Figure [8} 


GEANT4 
FLUKA rescaled 


9.6 


9.7x10'' 


3.2xl0~* 


1.2x10^ 
«1.2xl0* 


4.0x10"** 
Ri4.0xl0-'^ 


^ 0.05 X Njco 



TABLE IV: Cosmogenic background as a function of assumed depth: estimated number of neutrons (with kinetic energy above 
20 MeV), neutral kaons and lambda's entering the detector per year, and normalised on the average number of muons entering 
the active volume per 10 ms, produced in cosmic muon interactions (hadrons accompanied by a detected muon inside LAr 
imaging have been vetoed). 
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FIG. 11: Energy distribution of lambda particles entering the 
LAr volume as predicted by FLUKA. 



erage distance of one meter from it. The possibility to 
veto cosmic muon passing in the vicinity of the detector 
has therefore been investigated. 

We considered annular planes of active muon detectors 
with an inner radius of 35 m and an outer radius of 38 m, 
to detect muons that pass within 3 m from the LAr ac- 
tive volume (also here we have neglected the details of 
the liquid Argon tank). Two configurations have been 
taken into account: (a) two planes at the top and at the 
bottom of the LAr detector; (b) three planes as shown 
in Figure [T] In the simulations no assumptions are made 



regarding their composition, they have a thickness of 1 
cm and 100% efficiency for muon crossing the counters. 
In practice, they could be realized with Resistive Plate 
Chambers (RPC) 0, [Ml or Glass RPCs (GSC) [9^. 
The total area needed is about 3800 m^ per plane. In 
the MONOLITH proposal [9^ the total active area was 
envisaged to be « 54000 m^. Therefore GSC had to 
be designed in order to allow large-scale production. For 
techniques on the industrial large area production of glass 
RPCs see e.g. ReL ||. 

From the simulation in the under the hill configura- 
tion, we obtained that about 60% of the muons enter the 
veto region (i.e. a volume 20 m high and 3 m wide out- 
side the LAr detector) from the top. This value was used 
as a normalisation to calculate the fraction of muons de- 
tected with the two veto configurations, and the results 
are stated in Table |Vl Even with the two veto planes 
configuration, more than 80% of the muons could be de- 
tected. In conclusion, an active very large area veto com- 
posed of 2 or 3 planes could suppress cosmogenic back- 
ground by an order of magnitude. In the flat Earth sur- 
face geometry or at larger depths, we expect the suppres- 
sion to be even stronger since muons are more vertically 
distributed than in the hill configuration. These results 
will be used in the next section to compute sensitivities 
to nucleon decays. 



IV. NUCLEON DECAY ANALYSIS 

The list of studied decay channels in summarized in 
Table |Vll 

The approach to discriminate between signal and back- 
ground is based on a set of sequential selection cuts. Final 
state topology and event kinematics provide the selection 
criteria. At first, we ignore the cosmic muon- induced 
background, and apply cuts until the atmospheric neu- 
trino background can be considered as irreducible. 
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FIG. 12: Radial distribution of photo-nuclear interaction ver- 
texes producing a neutron that enters the detector at a depth 
of 3 km w.e. 



FIG. 13: Scatter plot showing the regions populated by fully 
simulated muons, kaons and protons on the < dE/dx > vs. 
kinetic energy plane. 



Configuration 


Muons (1.6x10®/ 
in 1"' veto plane 


I simulated per 
in veto volume 


configuration) 
% detected 


Full volume 


228 


377 


100% 


3 veto planes 


231 


361 


94% 


2 veto planes 


243 


343 


85% 



TABLE V: Fraction of muons entering the veto volume and 
detected using veto planes in different configurations. The 
veto volume is defined as a cylinder 20 m high and 3 meter 
wider in radius than the LAr volume. 



An efficient background reduction demands a good 
particle identification. Results based on a dedicated 
analysis show that the tagging efficiency of pions, kaons 
and protons is above 99% with contamination below 1%. 
In addition, the muon-pion misidentification is around 
40%. The analysis, based on fully-simulated Monte-Carlo 
events with single particles, combines several variables in 
a Multi-layer perceptron Neural Network architecture: 
the fit to the particle hypothesis on the < dE/dx > 
vs kinetic energy plane, the fitted particle mass and the 
energy released after the particle decay. Figure [T^ shows 
how, thanks to the fine detector granularity, particles of 
different species get clearly separated. Electron and pho- 
ton identification is based on an algorithm that distin- 
guishes single m.i.p. signals versus double m.i.p. signals 
by using the first hits of each identified electromagnetic 
track. 

The list of selection cuts to reduce the atmospheric 
neutrino contamination, final efficiencies and background 



estimations are described in sections IIV Al and IIVBI 

We then evaluate the backgrounds due to cosmogenic 
events. Our strategy is to reduce the contribution from 
this background at the level of the irreducible atmo- 
spheric neutrino background. The amount of cosmogenic 
background is strongly correlated with the travel path in- 
side the detector. Thus we will cut on the detector fidu- 
cial volume, while keeping the same topology and kine- 
matics criteria used for the neutrino-induced background 
(in this way, signal selection efficiencies are unaltered). 
The details on the cosmic muon-induced background es- 
timation are given in Section [IV CI 



A. Proton decay channels 

The sequential cuts applied for each channel are briefly 
described in the following paragraphs. The detailed list of 
cuts for the considered proton decay channels are listed 
in Table IVIII Survival fraction of signal (first column) 
and background events from the different atmospheric 
neutrino interactions after selection cuts are applied in 
succession are also listed. Backgrounds are normalized to 
an exposure of 1 Mtonxyear. The final efficiencies and 
expected background events after cuts are reported in Ta- 
ble IVII The published efficiencies, backgrounds and re- 
sults obtained by Super-Kamiokande are shown for com- 
parison. 

• p ^ e+ tt'^ channel: A simulated decay event is 
shown in Figure 1141 The two photons (from the 
tt'^ decay) and the positron flying in the opposite 



15 









This paper (LAr TPC) 


Super-Kamiokande resuits [5, 101 






Decay 
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mode 


(%) 


bacl:ground 


(%) 


background 
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100 ktonxyear 




92 ktonxyear 


90% C.L. 


(pl)p 






TT 


45.3 


0.1 


40 


0.2 


1.6 xlO^^ 


(p2)p 




U 


41.9 


78.2 








(p3)p 




V 


96.8 


0.1 
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6.0 (A'+ ^ TV+TV") 


0.7 
0.6 


2.3 xlO^'^ 


(p4) p 






44.8 


0.8 


32 
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(p5) p 
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0.04 


1.0 xlO^^ 
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7 


98.0 
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73 
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7 


98.0 


< 0.2 


51 


0.2 




(P9) P 




7r+ 7f+ 


97.6 


0.1 








(plO) p^e 




18.6 


2.5 








(ill) n 




V 


45.1 


47.4 








(n2) n 




K+ 


96.0 


< 0.2 








(n3) n 






44.4 


0.8 








(n4) n 




7V + 


44.8 
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TABLE VI: Summary of studied decay modes: signal detection efficiency and expected atmospheric neutrino background 
(normalized to 100 ktonxyear exposure) after selection cuts. Where available, the efficiencies and background results of Super- 
Kamiokande are given for comparison. The published results obtained by Super-Kamiokande are shown for completeness [sl-fioj]. 



direction are clearly visible. In the chosen readout 
view, the event spreads over about 120x100 cm^. 
Figure [12] shows the distributions of the follow- 
ing reconstructed kinematical quantities: the elec- 
tron momentum, the total momentum imbalance, 
the invariant mass and the total energy. The dis- 
tributions are split into the case where the pion 
leaves the nucleus (full histograms) and the case 
where it is absorbed (empty histograms). The elec- 
tron momentum histrogram has an arrow placed at 
460 MeV/c, the expected value without Fermi mo- 
tion and detector effects. It seems clear that two 
different set of cuts could be implemented to opti- 
mize the signal over background ratio in both cases. 
However, an attempt to look for "inclusive" decays 
p e+(7r°) without condition on the 7r°, yields 
order of magnitudes worse background conditions, 
and was not considered further. 

The list of cuts is presented in Table IVIII The 
idea is to have a balanced event, with all particles 
identified as such, and with a total visible energy 
close to the proton mass (see Figure [TCI). Only 
one background event for 1 Mtonxyear exposure 
survives the cuts, for a signal efficiency of about 
45%. 

• p — * 7r+ channel: Almost 45% of the events that 



belong to this channel can not be detected since 
the 7r+ gets absorbed by the nucleus. The cuts 
are based on the requirement of absence of charged 
leptons, protons, neutral pions, the presence of one 
charged pion and a total energy between 350 and 
650 MeV. The result is that, for a ~42% efficiency, 
the expected background at 1 ktonxyear exposure 
is '^0.8 events. 

• p — > 9 channel: this is a quite clean chan- 
nel due to the presence of a strange meson and no 
other particle in the final state (see Figure [T7)) . The 
kaon particle identification is performed and apply- 
ing the cuts listed in Table IVIII yields an efficiency 
'^97% for a negligible background. The correlation 
of the reconstructed invariant mass and total mo- 
mentum is shown in Figure [TSl 

• p ^ 1.1^ tt" channel: Almost 53% of the times 
the tt'^ is detected. In this case, cuts are similar 
to the e+Tr" channel. The efficiency remains high 
('^45%), while the background is ~ 8 events for 
1 Mtonxyear. 

• p — > e+ and p — > /i+ channels: We con- 
centrate on final states having a Kg, since a large 
fraction of the will leave the detector without 
decaying or will suffer hadronic interactions. In 
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Efficiency (%) 


Atmospheric neutrino 


sources 


Cuts 


(pi) p^e+ TT° 


ly^ CC 


Ue CC 


CC 


CC 


u NC 


iy NC 


One 7r° 


54.0% 


6604 


2135 


15259 


5794 


8095 


3103 


One e-shower + no other charged tracks 


50.9% 


1188 


656 


1 











ptot < 0.4 GeV 


46.7% 


454 


127 














0.86 GeV < E„i, < 0.95 GeV 


45.3% 


1 

















Cuts 


(p2) p ^ t:'^ P 


CC 


CC 


/^^ CC 


CC 


u NC 


ly NC 


No e-shower, no muons, no 7r° 


92.6% 








34 





56515 


26482 


One charged pion 


55.7% 








8 





5632 


2027 


No protons 


50.0% 








4 





2930 


1136 


0.35 GeV < Total E < 0.65 GeV 


41.9% 








2 





605 


175 


Cuts 


(p3) p^K+ i? 


i/e CC 


l^e CC 


I'm CC 


I^M CC 


ly NC 


ly NC 


One kaon 


96.8% 


308 


36 


871 


146 


282 


77 


No other charged tracks, no tt" 


96.8% 














57 


9 


E„i, < 0.8 GeV 


96.8% 














1 





Cuts 


(p4) p II -" 


CC 


;7, CC 


;/„ CC 


fy,, CC 


;/ NC 


fy NC 


One muon, one 


52.8% 








11334 


4452 








No protons, no charged pions 


50.0% 








1754 


1369 








0.86 GeV < Total E < 0.93 GeV 


45.3% 








64 


41 








Total Momentum < 0.5 GeV 


44.8% 








5 


3 








Cuts 


(p5) p^^i+ K°s 


J/e CC 


Ue CC 


CC 


CC 


u NC 


iy NC 


One muon + 2 charged or 2 neutral pions 


100% 


8178 


2771 


106861 


27274 


7099 


2540 


0.4 < Invariant mass of pions < 0.6 GcV 


97% 








5 


8 


6 


2 


Pioi < 0.() Gi'\' 


9:.i.l% 




















Cuts 


(p6) p^e+K° 


CC 


CC 


CC 


I^M CC 


ly NC 


iy NC 


One e-shower + 2 charged or 2 neutral pions 


100% 


59759 


11673 


31 





2 


1 


0.4 < Invariant mass of pions < 0.6 GeV 


97.0% 


2 


2 














Ptot < 0.6 GeV 


94.0% 




















Cuts 


(p7) p ^ e+ 7 


Ue CC 


i^e CC 


CC 


CC 


ly NC 


ly NC 


One e-shower + no other charged track 


100% 


32434 


6837 














Only one photon 


99.0% 


110 


11 














> 0.2 GeV 


98.0% 




















Cuts 


(p8) p ^ /i+ 7 


CC 


CC 


I'm CC 


CC 


ly NC 


ly NC 


One muon -|- no other charged track 


100% 


5302 


1878 


54889 


15872 


4680 


1764 


Only one photon 


99.0% 


7 


4 


164 


13 


9 


7 


p^ > 0.2 GeV 


98.0% 




















Cuts 


(p9) p ^ /x" 7r+ K+ 


CC 


l^e CC 


/^M CC 


CC 


ly NC 


[y NC 


One Kaon 


98.8% 


308 


36 


871 


146 


282 


77 


One muon 


98.2% 


1 





867 


146 








No e- showers 


98.2% 








844 


145 








0.6 GeV < Total E < 1 GeV 


97.6% 








1 











Cuts 


(plO) p — > e"*" Tr"*" 7r~ 


CC 


CC 


I'm CC 


CC 


ly NC 


ly NC 


One e-shower, no muons 


100% 


59755 


11673 














Two charged pions, no protons 


19.4% 


714 


302 














0.65 GeV < Total E < 1 GeV 


19.0% 


33 


8 














Total Momentum < 0.57 GeV 


18.6% 


21 


4 















TABLE VII: Detailed list of cuts for the considered proton decay channels. Survival fraction of signal (first column) and 
background events through event selections applied in succession. Backgrounds are normalized to an exposure of 1 Mtonxyear. 
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FIG. 14: Simulated p e n event. The displayed area 
covers 120 x 100 cm^. 



FIG. 16: Kinematic cut in the p e n channel: in the 
plane defined by the invariant mass and the total momen- 
tum, crosses represent background and spots signal events. 
The band indicates the cut region (0.86 GeV < Total En- 
ergy < 0.95 GeV), i.e. all events inside the band are accepted. 



p — ^ Tc" Kinematics 




0.2 0.4 0.6 O.S 1 0.4 0.6 0.8 1 



Invariant Mass (GeV) Measured Energy (GeV) ^IG. 17: Simulated p ^ K+ u event. The displayed area 

covers 34 x 90 cm'^. 



FIG. 15: p — > e n channel: Distributions of some kinematic 
variables for the exclusive (full histograms) and the inclusive 
(empty histograms) scenarios. The arrow in the first plot 
indicates the value that would have the positron momentum 
if no Fermi motion and no detector effects were present. 



Figure [13 simulated events of p ^ e+ + and 
p — > /i+ + K'g decays are shown, where the Kg 
decays into charged pions. A p — > /Lt+ -I- de- 
cay where K" decays into neutral pions is shown 
in Figure [201 A sophisticated treatment to re- 
cover K^ events has been neglected at this stage. 
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(0.65<-B^is<l GeV) complemented with a cut on 
the total momentum {ptot<570 MeV/c) are suffi- 
cient to reduce the contamination of 25 events for 
a 1 Mton xyear exposure. On the other hand, the 
cuts remove less than 1% of the signal events. 



B. Neutron decay channels 

The sequential cuts applied for each channel are briefly 
described in the following paragraphs. The detailed list of 
cuts for the considered neutron decay channels are listed 
in Table IVlTll Survival fraction of signal (first column) 
and background events from the different atmospheric 
neutrino interactions after selection cuts are applied in 
succession are also listed. Backgrounds are normalized 
to an exposure of 1 Mton xyear. The final efficiencies 
and expected background events after cuts are reported 
in Table IVD 



FIG. 18: Kinematic cut in the p — > K iy channel: in the 
plane defined by the invariant mass and the total momentum, 
crosses represent background and spots signal events. The 
band indicates the cut region (Total Energy < 0.8 GeV), i.e. 
all events in the band are accepted. 



K° 's mainly decay to two pions (either charged or 
neutral). Simple requirements (an identified lep- 
ton in the final state accompanied by two charged 
(neutral) pions, invariant mass of the pion system 
consistent with the K" mass and total momentum 
below 0.6 GeV) result in a negligible background 
contamination for signal selection efficiencies above 
90% (see Tabic IVUl. The ovcraU efficiency for the 
p n+ie+) K° channel is therefore 46.7% (47%). 

p e^^ and p — » /i^7 channels: These channels 
provide a very clean signal thanks to efficient elec- 
tron and photon separation (see e.g. Ref. [30] for a 
discussion in the context of e/Tr" separation). The 
simple final event topology (a single charged lep- 
ton accompanied by an energetic photon) allows 
to reduce the expected background to a negligi- 
ble level while keeping a signal selection efficiency 
close to 100% (see Table IVlH) . Simulated events of 



p ^ e'^ + ^ decay and p 
Figure im 



7 are shown in 



• p — > /i^ 7r+ channel: by tagging the presence 
of a and a p" , the background is reduced at 
the level of ~1 event for 1 Mtonxyear exposure. 
A final cut on the total visible energy removes any 
background for a signal of ~ 97%. 

• p ^ 7r+ TT^ channel: The most favorable 
scenario occurs when the three particles are de- 
tected. A tight cut on the total visible energy 



• n ^ -K^ V channel: About ~45% of the signal 
events are "invisible" because the 7r° is absorbed 
in the nucleus. The rest can be disentangled from 
the background by cutting on the total visible en- 
ergy and on the total momentum (>0.35 GeV). The 
final efficiency is ~45% for ~0.5 background events 
at 1 kton xyear exposure. 

• n ^ e~ channel: We profit here from the pres- 
ence of one kaon, one electron and the absence of 
muons and pions in the final state. These require- 
ments, together with a loose cut on the visible en- 
ergy (0.75<ii^«is<0.95 GeV), eliminates the back- 
ground for an almost untouched efficiency. 

• n — > e+ vr^ channel: This channel is similar to 
the p ^ tt" previously reported. When the 7r~ 
is detected, a good efficiency can be reached for a 
negligible background by applying two simple cuts 
to bound the total energy {0.75<Evis<l GeV) and 
the positron momentum (0.35<pe<0.6 GeV). 

• n — > fi^ 7r+ channel: This channel is treated in a 
similar way to the n ^ e+ tt^. The distribution 
of the total momentum for signal and background 
events is shown in Figure [22] and the cuts in Ta- 
ble IVIIII We require the presence of a pion on the 
final state and a cut on the total energy. The plot 
shows the position of the last cut on total momen- 
tum. The final efficiency is ^^45% for the exclusive 
channel and almost no background events expected. 



C. Cosmic muon-induced background estimation 



As described in Section IIII Dl the cosmogenic back- 
ground has been computed in three steps. In the third 
step, neutrons, neutral kaons and lambda's entering the 
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FIG. 19: Simulated p ^ + (left) and p — + fi'^ + K'^ (right) decays. The neutral kaons decay into two charged pions. 




detector were generated according to their expected en- 
ergy spectrum. We simulated 2x10^^ events for each par- 
ticle species and used values in Table HVl for normaliza- 
tion. The neutral particles are propagated inside the liq- 
uid Argon volume until they interact inelastically or de- 



cay. The secondary particles produced in those processes 
were used in the nucleon decay analysis. 

As concrete examples, we discuss in detail the two pro- 
ton decay channels, p D and p tt"*" 9, and one 
neutron decay channel n — > tt'^ j/. Other decay chan- 
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Cuts 


Efficiency (%) 


Atmospheric neutrino sources 


(nl) n — > tt'' 1/ 


CC 


CC 


ly^ CC 


CC 


u NC 


jy NC 


One 7r° 


56.2% 


6604 


2135 


15259 


5794 


8095 


3103 


No muons, no electrons, no charged pions 


56.1% 








2 





4722 


1840 


No protons 


52.6% 














2964 


1184 


0.35 GeV < Total E < 0.55 GeV 


45.4% 














469 


181 


Total Momentum > 0.35 GeV 


45.1% 














362 


112 


Cuts 


(n2) n^e- K+ 


Ue CC 


i^e CC 


CC 


u^. CC 


u NC 


jy NC 


One e-shower, one kaon 


97.0% 


299 


36 


11 











No TT° , no muons 


97.0% 


138 


14 














No charged pions 


97.0% 


80 


5 














0.75 GeV < Total E < 0.95 GeV 


96.0% 




















Cuts 


(n3) n e'^ 


CC 


CC 


/^M CC 


CC 


u NC 


iy NC 


One e-shower, one charged pion 


59.6% 


8137 


2755 


6 











No tt'^ , no muons, no protons 


57.4% 


3855 


1282 














0.75 GeV < Total E < 1 GeV 


52.4% 


499 


187 














0.35 GeV < Ppositron < 0.6 GeV 


51.3% 


216 


73 














Total Momentum < 0.35 GeV 


44.4% 


7 


1 














Cuts 


(n4) n ^ jj.^ 7r+ 


CC 


CC 


CC 


CC 


u NC 


iy NC 


One muon, one charged pion 


59.4% 


1559 


454 


15931 


6569 


2291 


1055 


No tt", no e-shower. No protons 


53.6% 








7830 


2924 


824 


444 


0.8 GeV < E„i, < 1.05 GeV 


49.8% 








1064 


408 


137 


56 


ptot < 0.35 GeV 


44.8% 








18 


2 


5 


1 



TABLE VIII: Detailed list of cuts for the considered neutron decay channels. Survival fraction of signal (first column) and 
background events through event selections applied in succession. Backgrounds are normalized to an exposure of 1 Mtonxyear. 



nels with a lepton in the final state or more constrained 
kinematics and topologies, will be less affected by cos- 
mogenic backgrounds. For instance, the probability to 
misidentify charged pions as muons is at these momenta 
typically less than 10~^. 

In the analysis, we started by applying loose kinemat- 
ical cuts that rejected cosmogenic events clearly incom- 
patible with nucleon decay signals [95]. For instance, for 
the channel p — > tt"*" the expected pion energy in an 
ideal detector would be 0.48 GeV. We accepted charged 
pions with energy in the range 0.35-0.65 GeV, since the 
measured energy is smeared by Fermi motion and de- 
tector effects. In the case of strange mesons, we also 
required an identified kaon in the final state and a to- 
tal energy below 0.8 GeV. In a second step, same cuts 
as those defined in Tables IVIII and IVIIII were applied to 
estimate final background contaminations. 

Table IIXI summarizes for each considered channel the 
remaining muon-induced background level after this se- 
lection. The contamination coming from neutrons, kaons 
and lambdas interactions at different detector depths are 
shown. We found that the tf^ and tt*' backgrounds are 
predominantly produced by neutrons, and the back- 
ground is dominated by the entering into the detec- 
tor. For some channels we obtain zero events after the 
final cut, hence an upper limit for these events occur- 



ring in the detector is given. The A induced background 
events appeared to be negligible in comparison to the 
other sources. 

The photo-nuclear interactions, from which at least one 
neutral particle enters the detector, occur mostly at an 
average distance of one meter from the detector walls (see 
Figure fT2)). Moreover, the number of background events 
reduces itself essentially exponentially along their path 
inside the detector due to the self shielding properties 
of Argon. Therefore, an obvious action to reduce back- 
ground consists on cutting on the sides of the detector, 
excluding events produced at a distance smaller than d 
from the wall (see Figure 1^5)) . The fiducial volume is re- 
duced accordingly. The value for the cut distance d for 
each detector depth can be chosen in such a way that the 
remaining muon-induced background is of the same order 
than the irreducible atmospheric neutrino background. 

In order to compute the effect of the distance cut on 
the background events, we consider here the most con- 
servative situation, i.e. the incident neutral background 
source particles {n and iC°) move horizontally from the 
top/bottom side wall into the detector (the shortest pos- 
sible way). After reaching the cut distance d, the parti- 
cles move along the diagonal of the inner detector cylin- 
der (the longest possible way). We obtain then a simpli- 
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Depth 


Channel 


Backgroi 
(part 
Neutron 


md sour 
icles/yef 
K° 


ce N'il 
ir) 

A 


Cosmogf 
Distance cut 
d (m) 


3nic backgroun 
Fiducial mass 
(kton) 


d reduction 
Background Nt, 
(events/year) 


~0.5 km w.e. 
(188 m rock) 
FLUKA 


p ■K^h' 

n -K^v 
p K+u 


570 
450 


135 


8 


1.5 
1.7 
6.6 


92 
91 
66 


76 
46 
0.1 


~ 1 km w.e. 
(377 m rock) 
FLUKA 


P TX^U 

p K+u 


200 
130 


39 


2.3 


0.7 
0.75 
5.45 


96 
96 
71 


77 
47 
0.1 


~ 3 km w.e. 
(1.13 km rock) 
FLUKA 


P TX^U 

n TT^'iy 
p K+u 


4.0 
2.6 


0.74 


- 





1.8 


100 
100 
90 


4.0 
2.6 
0.1 


Under the hill 
(see Figure |8j 
GEANT4 


P TX^V 

n TT^'iy 
p -> K+u 


2900 
2300 


36-360 


- 


2.7 
2.9 
5.4-7.5 


85 
84 
72-62 


76 
46 
0.1 


Under the hill 
+ two veto planes 
GEANT4 


P TX^V 

n Tx^u 
p K+v 


430 
340 


5-54 


- 


1.3 
1.5 
3.65-5.75 


93 
92 
80-70 


76 
46 
0.1 


Under the hill 

+ three veto planes 

GEANT4 


P —t TX^V 

n — > Tx^i? 
p K+v 


170 
140 


2-20 




0.6 
0.8 
2.8-5 


97 
95 

85-74 


77 
46 
0.1 



TABLE IX: Cosmogenic background for three selected channels: estimated number of background events per year that survive 
a kinematic selection. The contamination coming from neutrons, kaons and lambdas interactions at different detector depths 
are shown. For each detector depth, the radial cut distance and the final fiducial volume to reduce cosmogenic background to 
the level of the irreducible atmospheric background (resp. 78.2 for p — > tx+v, 47.4 for n — > tx^v and 0.1 for p — > K+v for an 
exposure of 100 kton x year) is listed. The range for kaon background is reflecting uncertainty on kaon yields due to differences 
between FLUKA and GEANT4 results. 



fied expression for the normalised background events: 



{R - d)2 



(4) 



where and Nj, are the numbers of events before and 
after the cut on the fiducial volume, and cx is an effective 
interaction length found to be ~ 71 cm for neutrons and 
^ 86 cm for kaons from the MC simulations. With the 
given number of neutrino produced background events 
as Nb, the cut distance d can be calculated from the 
following equation: 



d^R + 2caW 




(5) 



where the W{x) is the Lambert W-function (or Omega 
function, inverse of f{w) — we"^). 

Table Hxl summarizes . for each considered channel, the 
distance cuts applied together with the remaining fidu- 
cial mass. For each detector depth, the radial cut dis- 
tance and the final fiducial volume to reduce cosmogenic 
background to the level of the irreducible atmospheric 
background (resp. 78.2 for p tt+v, 47.4 for n — > 7r°P 
and 0.1 for p K+v for an exposure of 100 ktonxyear) 



is listed. The range for kaon background is reflecting un- 
certainty on kaon yields due to differences between FLUKA 
and GEANT4 results. 

We observe that the background due to charged and 
neutral pions can be reduced to the same level of expected 
neutrino-induced backgrounds without a big loss of the 
detector fiducial mass, even at shallow depth (1 km w.e. 
or equivalently 380 m of ordinary rock). The rejection 
of strange mesons requires harder cuts. For a 1 km w.e. 
depth, a reduction of 30% of the mass is expected. How- 
ever the prospects for a shallow depth experiment are 
promising and good sensitivities for most of the analyzed 
channels are expected for an underground location of 1 
km w.e. (see Table IX)) . Similar conclusions hold when a 
0.5 km w.e. depth (about 200 m of ordinary rock) is con- 
sidered. Clearly the situation worsens very much if we 
consider a position close to surface. For a 50 m rock over- 
burden, the reduction of the pion background requires a 
reduction by more than 30% of the usable Argon. The 
reduction of the kaon background implies that half of the 
available Argon should be used as shield. 

Note that using an annular veto system to detect 
muons as explained in Section IIII El the induced back- 
ground can be reduced by a factor between ~ 5 and 
15 (according to the number of planes used), there- 
fore a looser cut on the distance d would be needed and 
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FIG. 21: Simulated p e'^ +'y and p ^ + 7 decays. 



it would be possible to have a bigger fiducial volume. For 
example, the under the hill configuration with three veto 
planes consents to reduce the distance cut from w 6 m 
to « 3 m (see Table E]). 



V. SENSITIVITY TO NUCLEON DECAY AND 
COMPARISON TO SUPER-KAMIOKANDE 
RESULTS 

We summarize here the the results from the previous 
sections for all channels (see Table rvT| and compute par- 
tial lifetime sensitivities. 

In case no signal is observed, limits to proton and neu- 
tron partial lifetimes t/B will be obtained using : 



ir/B)proton > ^ X T X e X 10^2 

3 3 

{T/B)neutron > — X T X e X 10^^ 



(years) (6) 



(years) (7) 



T is the exposure in kilotons x year, e the signal selection 
efficiency, and S the constrained 90% CL upper limit on 
the number of observed signal events, taking into account 
the fact that there are about 2.7 xlO'^^ protons and 3.3 
X 10^2 neutrons in 1 kton of Argon. S is found by solving 
the equation [ssj : 



E7=oPin,b + S) 



> 

s 



selected region 



I I Vj,CC + v^CC Background 
I I n — > |i" tC Signal 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Total Momentum (GeV) 



FIG. 22: n ^ jj. n channel: distribution of total momen- 
tum for events surviving the sixth cut of Table IVlnl 




(8) 



FIG. 23: Schematic view of the radial distance cut to suppress 
cosmogenic neutral particle background entering the sides of 
the detector. 



P(n,/i) is the Poisson function, b the estimated back- 
ground, a = 0.1 for a 90% CL, and, uq is equal to the 
closest integer number to b when computing the "detec- 
tor sensitivity" . 

For each nucleon decay channel we have computed the 
{t/B) limits as a function of the exposure. This is done 
by rescaling the number of expected background events 
at each exposure, and computing the corresponding up- 
per hmit (5*). The result, together with the detection 
signal efficiency (e), the total expected background and 
the final detector mass (after fiducial cuts) are reported 
in Table 1x1 In the channels p — > tt+i/, p K~^D, and 
n 7r°P, the results at different detector depths are 
given separately. 

The sensitivity of (t/B) for protons and neutrons as a 



23 



Channel 


Cut 


Total background 


r/B limit (years) 


t/B limit (years) 










efficiency (%) 


per year in fiducial volume 


1 year exposure 


10 years exposure 


(pl)p 




e+ 




45.3 


0.1 


0.5 


xlO^* 


0.4 


xlO^^ 


(P2)p 




7r+ 


V 


41.9 


82 (3 km w.e.) 


0.7 


xlO^^ 


0.3 


xlO^^ 












151 (1 km w.e.) 


0.5 


xlO^^ 


0.2 


xl034 












148 (0.5 km w.e.) 


0.5 


xlO^^ 


0.2 


xlO^^ 












143 (Under the hill) 


0.4 


xlO^^ 


0.2 


xlO^^ 












149 (Under the hill+2 veto planes) 


0.5 


xlO^^ 


0.2 


xlO^"* 












152 (Under the hill+3 veto planes) 


0.5 


xlO^^ 


0.2 


xlO^^ 


(p3)p 




K+ V 


96.8 


0.2 (3 km w.e.) 


1.0 


xlO^^ 


0.6 


xlO^' 












0.2 (1 km w.e.) 


0.8 


xl034 


0.6 


xlO-^'' 












0.2 (0.5 km w.e.) 


0.8 


xlO^* 


0.4 


xlO^^ 












0.2 (Under the hill) 


(0.8-0.7) xlO'''' 


(0.5-0.4) xlO^^ 












0.2 (Under the hill+ 2 veto planes) 


(0.9-0.8) xlO'''' 


(0.5-0.5) xlO^^ 












0.2 (Under the hill+ 3 veto planes) 


(1.0-0 


8) xlO^"* 


(0.6-0.5) xlO^^ 


(p4) p 




^+ 


ttO 


44.8 


0.8 


0.4 


xl034 


0.2 


xlO^^ 


(p5) p 






7f» 


46.7 


< 0.2 


0.5 


xlO^'' 


0.5 


xlO^^ 


(p6) p 




e+ 


if" 


47.0 


< 0.2 


0.5 


xl034 


0.5 


xlO^^ 


(p7)p 




e+ 


7 


98.0 


< 0.2 


1.1 


xl0^4 


1.1 


xlO^^ 


(p8)p 






7 


98.0 


< 0.2 


1.1 


xl034 


1.1 


xlO^^ 


(p9)p 








97.6 


0.1 


1.1 


xlO^* 


0.8 


xlO^^ 


(plO) 1 






18.6 


2.5 


0.1 


xl034 


0.5 


X lO'^"* 


(nl) n 






V 


45.1 


50 (3 km w.e.) 


0.1 


xlO^* 


0.5 


xl0^4 












92 (1 km w.e.) 


0.1 


xlO^'-' 


0.4 


xl034 












89 (0.5 km w.e.) 


0.1 


xlO^'' 


0.4 


xlO^* 












86 (Under the hill) 


0.1 


xlO^'' 


0.3 


xl034 












90 (Under the hill+ 2 veto planes) 


0.1 


xl034 


0.4 


xlO^"* 












91 (Under the hill+ 3 veto planes) 


0.1 


xlO^'' 


0.4 


xlO^"* 


(n2) n 






K+ 


96.0 


< 0.2 


1.4 


xlO^* 


1.4 


xlO^^ 


(n3) n 


— ^ 


e+ 




44.4 


0.8 


0.4 


xl034 


0.2 


xlO^^ 


(n4) n 






7r+ 


44.8 


2.6 


0.4 


xlO=''' 


0.2 


xlO^^ 



TABLE X: Summary table of the main nucleon decay search results. For each channel, the signal detection efficiency and the 
total expected background are given together with the lifetime limit at 90% CL for 1 and 10 years of data-taking. For the 
p K^u the range reflects the difference of yields of kaons per muon from FLUKA and GEANT4. 



function of the exposure is illustrated in Figure [24] con- 
sidering only atmospheric neutrino background (left) and 
including the cosmogenic background (right). In the lat- 
ter case the plot is shown as a function of exposure in 
years since there is a reduction in fiducial mass taken 
into account to reduce the background according to Ta- 
ble IIXI while the sensitivity curves without cosmogenic 
background are plotted as a function of the exposure in 
kton X year. 

The liquid Argon TPC, offering good granularity and 
energy resolution, low particle detection threshold, and 
excellent background discrimination, can yield good sig- 
nal over background ratios in many possible decay modes, 
allowing to reach partial lifetime sensitivities in the range 
of 10'*'* — 10'^^ years often in background free conditions 
up to exposures of 1000 kton x year. This situation is 
optimal for discoveries at the few-events level. In partic- 



ular: 

• Multi-prong decay modes like e.g. p —>■ ^i^n^K^ 
or p — > e~^iT'^TT~ and channels involving kaons like 
e.g. p — > K~^i?, p e'^K'^ and p — > ^'^K^ are par- 
ticularly suitable, since liquid Argon imaging pro- 
vides typically an order of magnitude improvement 
in efhciencies for similar or better background con- 
ditions compared to Super-Kamiokande results. 

• Up to a factor 2 improvement in efhciency is ex- 
pected for modes like p — s- e+7 and p — > /i"''7 thanks 
to the clean photon identification and separation 
from TT*^. 

• Channels like p — s- e^Tr'^ or p — > /i^Tr*^, dominated 
by intrinsic nuclear effects, yield similar efficiencies 
and backgrounds as in Super-Kamiokande. 
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FIG. 24: Running of the nucleon partial lifetime sensitivity {t/B at 90% C.L.) as a function of exposure (left) considering only 
atmospheric neutrino background (right) with cosmogenic background and corresponding reduction in fiducial mass taken into 
account (see text). 



Thanks to the self-shielding and 3D-iniaging properties 
of the liquid Argon TPC, these results remains valid even 
at shallow depths where cosmogenic background sources 
are important. A very large area annular active muon 
veto shield could be used in order to further suppress 
cosmogenic backgrounds at shallow depths. For exam- 
ple, our results show that a three plane active veto at 
a shallow depth of about 200 m rock overburden in the 
under the hill configuration yields similar sensitivity for 
p — > K'^D as a 3 km w.e. deep detector. 

VI. CONCLUSIONS 

The most direct sign for Grand Unification is the ex- 
perimental detection of proton or bound-neutron decays. 
In order to reach partial lifetime in the relevant range, 
new generation massive underground detectors with fine 
tracking and excellent calorimetry are needed to suppress 
backgrounds with a good signal selection efficiency. Fur- 
thermore, the detector should be sensitive to several dif- 
ferent channels in order to better understand the nucleon 
decay mechanism. 

In this paper we have estimated the discovery poten- 
tial of liquid Argon TPC detectors and compared it to 
the existing results of Super-Kamiokande. We have ana- 
lyzed many possible neutron and proton decay channels 
and for the first time the two sources of background con- 
tamination: atmospheric neutrinos and cosmogenic back- 
grounds. 

Our results indicate that a liquid Argon TPC, if scaled 
at the relevant mass scale of the order of 100 kton, offer- 
ing good granularity and energy resolution, low particle 



detection threshold, and excellent background discrim- 
ination, can yield good signal over background ratios 
in many possible decay modes, allowing to reach par- 
tial lifetime sensitivities in the range of 10'^^ — 10^^ years 
with exposures up to 1000 ktonxyear, often in quasi- 
background-free conditions optimal for discoveries at the 
few events level, corresponding to atmospheric neutrino 
background rejections of the order of 10^. 

Multi-prong decay modes like e.g. p fi^n^K^ 
or p —> e^'K^TT~ and channels involving kaons like e.g. 
p — i- K^v, p — + e'^K^ and p — > are particularly 

suitable, since liquid Argon imaging provides typically an 
order of magnitude improvement in efficiencies for simi- 
lar or better background conditions compared to Super- 
Kamiokande. 

Up to a factor 2 improvement in efficiency is expected 
for modes like p — > e^7 and p — » /i^7 thanks to the clean 
photon identification and separation from tt*'. 

Channels like p e^Tr*^ or p — > /^^tt", dominated 
by intrinsic nuclear effects, yield similar efficiencies and 
backgrounds as in Super-Kamiokande. 

Thanks to the self-shielding and 3-D imaging proper- 
ties of the liquid Argon TPC, this result remains valid 
even at shallow depths where cosmogenic background 
sources are important. We consider the possibility of a 
very large area annular active muon veto shield in order 
to further suppress cosmogenic backgrounds at shallow 
depths. 

In conclusion, we find that this class of detectors does 
not necessarily require very deep underground labora- 
tories, like those typically encountered in existing or 
planned sites, to perform very sensitive nucleon decay 
searches. As a concrete example, we considered a three- 
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plane annular active veto surrounding a detector at a 
shallow depth of about 200 m rock overburden in an un- 
der the hill configuration and found that it yields similar 
sensitivities, in particular for p — K'^v, as for a 3 km w.e. 
deep detector. 

In addition to a successful completion of the required 
technological R&D necessary to reach the relevant mass 
scale of 100 kton in a cost-effective way, we point out the 
importance of an experimental verification of the liquid 
Argon TPC physics potentialities to detect, reconstruct 
and classify events in the relevant energy range. This ex- 
perimental verification will require in addition to possible 
specific tests in charged particle beams, the collection and 
study of neutrino events in GeV range with statistics of 



the order of lOO'OOO events or more. 
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The scintillation light signal has not being used in the 
present analysis, but will in the case of the actual exper- 
iment yield the To of the event, allowing to reconstruct 
the position of the event along the drift direction. 
The DIS process includes single-pion production via 
baryon resonance excitation within a parton model for- 
malism, inspired from the concept of duality. 
If the rock is well known and of very good quality in a 
given site, a new hole in virgin ground can be considered. 
See HH and references therein. 

The geography of the chosen location was estimated with 
the Google Earth tool. See ,http: //earth .google . com[ 



